Cumulative identical spin rotation effects in collisionless trapped atomic gases 
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We discuss the strong spin segregation in a dilute trapped Fermi gas recently observed by Du et 
al. with "anomalous" large time scale and amplitude. In a collisionless regime, the atoms oscillate 
rapidly in the trap and average the inhomogeneous external field in an energy dependent way, which 
controls their transverse spin precession frequency. During interactions between atoms with different 
spin directions, the identical spin rotation effect (ISRE) transfers atoms to the up or down spin state, 
depending on their motional energy. Since low energy atoms are closer to the center of the trap 
than high energy atoms, the final outcome is a strong correlation between spins and positions. 

PACS numbers: 



Spin waves in dilute gases were predicted at the be- 
of the eighties [H, [3] and confirmed experimen- 
tally soon after [3;, Jj, ^ . They can be understood in two 
equivalent ways, either as a consequence of spin mean 
field or in more microscopic terms as the cumula- 
tive result of the identical spin rotation effect (ISRE) - 
an effect taking place during binary collisions between 
identical atoms Q]. Similarly, the Faraday effect is a ro- 
tation of the spin of photons that can be seen, either as 
a consequence of a macroscopic index of refraction, or as 
resulting from the accumulation of microscopic forward 
scattering events between photons and atoms. 

Experiments with ultracold atomic gases have renewed 
the interest in spin waves. In 2002, a group at JILA 
showed that, in an trapped ultracold atomic gas of 
bosons, the ISRE can result in a spontaneous spatial 
"segregation" of two atomic internal states |1) and |2) 
(equivalent to a pseudo-spin 1/2). Several groups then 
proposed a theoretical explanation of these observations, 
using either one-dimensional spin 1/2 hydrodynamic [7] 
or kinetic [1, Q equations. More recently, Du et al. from 
Duke University Yu\ did an experiment that explores the 
properties of spin waves in quantum gases of fermions 
(^Li) in the collisionless Knudsen regime, while most pre- 
vious experiments were per formed in the hydrodynamic 
regime (see nevertheless [il|). The spin segregation they 
observe is a hundred times larger and a hundred times 
slower than would be predicted by hydrodynamic theory. 
They call this spectacular effect "anomalous spin segre- 
gation" and suggest that its explanation may require "a 
modification of spin wave theory or possibly a new mech- 
anism" for fermions. 

The purpose of this letter is twofold. First we argue 
that no modification of spin wave theory is necessary to 
understand the experiment; the physical mechanism be- 
hind the observations is the usual ISRE. The difference 
between bosons and fermions is not essential; what is 
important is the collisionless regime. Second, we discuss 
why this collisionless regime, combined with the presence 
of a trap, gives access to unexpected and interesting new 
physics. In fact, one can observe the ISRE almost as in 
an ideal experiment, where a single spin polarized atom 



is sent through a target of a gas polarized in another 
direction, and where the spin direction of the outgoing 
atom is measured. Moreover, when the trap potential 
sends back the atom through the same target, the rota- 
tions of the spins are additive (the symmetry of ISRE 
is the same as that of the Faraday effect). As a result, 
one reaches situations where the spins of the atoms are 
correlated to their motional energies in the trap, instead 
of their positions; these situations are inaccessible in the 
hydrodynamic regime, where the value of the spin cur- 
rent is determined locally by a thermal average over many 
collisions. 

We begin with a comparison between the JILA and 
Duke experiments. In the latter, the peak density n(0) 
and the typical scattering length ai2 are respectively 15 
and 20 times smaller than in the JILA experiment, mak- 
ing the diluteness factor l/rt(0)af2 of the gas 10^ times 
larger. The differences can be expressed in terms of 5 
relevant time scales: (1) the radial trap period 2iT/ujrad', 
(2) the axial trap period 27r/w; (3) the average time be- 
tween (lateral) collisions r ~ {A7ral2n{0)vT)~^ , where 
Vt = y^ksT/m is the thermal velocity, T the tempera- 
ture and TO the atomic mass; (4) the typical spin preces- 
sion period in the external magnetic field 27r/|r2o|; and 
(5) the typical precession period in the spin mean-field 
Tfwd ^ m/ h\ai2\n,{0) . These time scales [in ascending 
order for the Duke exp.] are compared in the table: 
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Duke exp. [1_0] 
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6.9 


100 


300 


5000 


JILA exp. [6] 


4 


143 


170 


14 


10 



where ai2 ^ — 5ao and rio/27r « —10 Hz are typical val- 
ues for the Duke exp. The ISRE is then strong: 
T/Tf^cL ~ h/ (mvT\a,i2\) ^ 1- In the JILA experiment, 
an hydrodynamic description was qualitatively correct 
because of the relatively small value of lot < 1. In the 
Duke experiment lot ~ 4500: the gas is so dilute that 
its dynamics is well described in the collisionless limit 
(lateral collisions are ignored). On average, an atom os- 
cillates 700 times in the axial trap between two collisions. 

We base our analysis on the formalism used in and 
fill, which gives a general frame for the study of non de- 
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generate quantum gases. It distinguishes two effects of 
the interactions: mean field effects (forward scattering in 
colhsions) and "real" collisions (lateral scattering), both 
with a full treatment of the effects of spin polarization 
and statistics. This leads to a kinetic equation with ar- 
bitrary position and momentum dependence. In the hy- 
drodynamic limit, the system remains close to local equi- 
librium, momenta can be integrated out, and the spin 
current is simply proportional to the local gradient of 
magnetization. But we do not make such an assumption 
and keep a full dependence on the variables. 

Consider a cigar-shaped harmonic trap with axial {x 
direction) frequency a;/27r and radial {y and z direc- 
tions) frequencies Wrad/STr. As w ^ tOradi an effec- 
tive ID description of the dynamics is possible. The 
peak density (per unit volume) n{0) is related to the to- 
tal number of atoms by = n{0){2n)'^^'^ xtt"^, where 
2 (resp. rj_ = ^/ ksT / mujl^j) is the 
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characteristic atomic cloud size in the axial (resp. radial) 
direction. The two internal states |1) and |2) are treated 
as an effective spin 1/2. The three relevant scattering 
lengths are ai2, an and 022; the corresponding coupling 
constants are Qij = ATrh?aij /m. In spin space, ey is the 
longitudinal unit vector, while e_Li, e±2 correspond to 
the transverse directions. The spin couples to an exter- 
nal effective magnetic field f2(x)e||, with an inhomogene- 
ity along x characterized by a curvature 17" (0). There- 
fore, the magnetic field is taken to be ^l{x) « Hqx'^/x'^, 
where fto = f2"(0)x^/2 is a characteristic spin precession 
frequency fl^. 

The kinetic equation for a spin 1/2 non-degenerate 
atomic gas (fermions or bosons) is written in terms of 
a 2 X 2 operator p obtained by a Wigner transform of the 
density operator: 

i \f{x,p,t)i + M{x,p,t)-a 



p{x,p,t) 



(1) 



Equivalently, we can reason in terms of the phase space 
density / and spin density M. In the Duke experiment, 
the gas is non-degenerate {T « 21 pK ^ Tp ~ 7pK). 
As in Q we make some approximations: both the den- 
sity mean-field |(7i2|n(0)/4 « /i x 2 Hz and the Stern 
and Gerlach segregation energies h\^lo\ « h x 10 Hz are 
negligible when compared to the harmonic confining po- 
tential miLj^x'^/2 = kBT/2 « /i x 280 kHz. Since lat- 
eral collisions are a small perturbation, we treat them in 
the simplest relaxation time approximation. Here a sin- 
gle coupling constant gi2 ^ appears, instead of three, 
because gu — (722 = results from the Pauli exclusion 
principle (collisions only occur between atoms in different 
spin states). We then obtain the equations Q: 

dtf = {dt + pd^ - xdp)f ~ ^LiIH (2) 



dtM. - [^{x)e\\ + egm] x M 



(3) 



where e = — 1 for fermions (and +1 for bosons); the su- 
perscript denotes local equilibrium phase space densi- 
ties; the dimensionless coupling constant g is defined as 



g = gi2n{0)/2Huj. Dimensionless units have been used: 
lengths are measured in units of xt, momenta in units 
of pt = mvT, angular frequencies in units of uj, times in 
units of 1/lu, and phase space densities / and M in units 
of n(0)/pT\/27r. The density n and the spin density m 
[both in units of n{0)] are defined as: 



n{x, t) = 



dp 



f and m(a;, t) = 



dp 
x/2^ 



M (4) 



such that 1 = / dxn/ y/2Tr [and e.g. e_Li — J dxm/y'ZTT 
if the gas is fully polarized in the e_Li direction] . The total 
density n = rii + n2 and the longitudinal spin density 
my = n2—ni are related to the internal state populations, 
while the transverse spin density m_L to their coherences. 

These equations are now solved numerically. Since the 
experiment starts with a 7r/2 pulse which suddenly trans- 
fers the spins to the transverse plane, the initial distri- 
butions are fo{x,p) = exp[— (x^ +p-^)/2] and Mo(a;,p) — 
fa{x,p)e±i. As /o solves the kinetic equation, the to- 
tal density does not evolve n{x) — exp(— x^/2) and the 
dynamics after the pulse can be expressed in terms of 
M only. The 3 dimensionless parameters appearing in 
the kinetic equations are taken as eg = 0.021 [corre- 
sponding to ai2 = — 4.6ao and e — —1], Qq — —0.07 
and 1/r = 2 x 10^'*. Fig. [T] compares the results 
with those plotted in fig. 2c of [l3|: both the maxi- 
mum segregation time {tmax ~ 300 ms) and amplitude 
(to|| (0, tma2;)/?^(0) ~ 30%) are in good agreement with 
the experimental data. Note also that the sign of the seg- 
regation is predicted correctly. It is clear that standard 
theory with a transport equation treated in the collision- 
less regime explains the observations without introduc- 
ing any new mechanism. The system eventually relaxes 
to equilibrium on the time scale of a few r, which can be 
very slow indeed because r ~ 5 s. Calculations for other 
values of ai2 are in fair agreement with the experimental 
data. In particular, changing the sign of ai2, of or of e 
simply reverses the role of the two internal states. When 
ai2 « 0, no significant segregation takes place. Fig. [2] 
shows the longitudinal spin density profile at saturation 
and a good agreement with the experiment results p^ . 

In order to better understand the numerical results of 
the previous section, approximate analytical solutions of 
the linearized [|.g| << 1] collisionless kinetic equations 
are now provided. In the absence of lateral collisions, the 
(axial) motional energy of an atom E = (x^ +p'^)/2 [in 
units of ksT] is a constant of the motion and M(x,p, t) 
conserves its norm. We therefore write My = /osin6', 
M_L = Mj_i+iMj_2 = fo cos 6*6"^, which define the trans- 
verse (j){x,p,t) and longitudinal 9{x,p,t) precession an- 
gles in phase space. Accordingly the initial conditions 
are 9{x,p,Q) — and 4'{x,p,Q) — 0. Eq. ([3]) translates 
into coupled equations of motion for 9 and 0: 



smf 



dt4> = r2 + egi{sm9) — 

cos 9 

dt9 ^ egluY{e">' {cos ee-">')} 



Re{e*'^(cos6'e-*'^)})(5) 
(6) 
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FIG. 1: Spin density at the center of the trap as a function of 
time m{x = 0, t) [in units of n(0)] for eg = 0.021, Oo = -0.07 
and 1/r — 2 x 10^*. We obtain a good agreement with the 
measurements of my = n2 — ni, taken from fig. 2c of [l3| and 
shown here with their error bars: both the amphtude and the 
time constants of the spin segregation are weU reproduced. 



verse spin precession angle 4>{x,p,t): 

^ = M\l+{t)x' + i7-(0/ - lo{t)xp\ (7) 



7o(t) = siT^t. The longitudinal 



where 7± (i) = t±- 
spin density m\\{x,t) is: 



£5 



dt' 



A Y Z Y 
— sm - arctan — — -r^-rr- 



(8) 



where Z[t') 



1 



7^(i')c)^, A{t,t') = 1 + (y)^, using the shorthand nota- 
tion s = sin {t — t') and c = cos (i — t'). 

Three regimes can be distinguished: (i) for times 
smaller than the trap period {t ^ 1), at a fixed posi- 
tion the transverse phase difference between two atoms 
with momentum p and p' is dominated by the x — p cor- 
relation term 70 (i) in eq. ([7]): 4>{x,p,t) — 4>{x,p\t) ~ 
not'^x{p — p') <C 1. This is the regime leading to correla- 
tion between velocity and transverse spin and giving: 



mn{x,t) w egnixfil - 2a;^)2f7o 




(9) 



FIG. 2: Numerical results (full line) for the longitudinal spin 
density at maximum segregation time m|| {x, tmax) [in units of 
n(0)] as a function of position; the parameters are the same 
as in fig. [T]and xt = 210 /xm. These results are in very good 
agreement with the measurements of [TO] (crosses) at the trap 
center and edges, without any adjustable parameter. The 
discrepancy close to x/xt ~ ±1 may be attributed to a small 
population imbalance in the experiment (~ 4% as estimated 
from the area under the experimental profile). The results of 
two approximations are also shown (dotted and dashed lines); 
see text for more details. 



where the momentum average is defined as (A) — 
J -^E=fQA for any function A{x,p,t). To first order in 

g, we can ignore the second term in eq. ^ and as- 
sume cos 9 ~ 1 in eq. ^ (furthermore ~ fo9 and 
M_L ~ /ne'*^). These linearized equations can be solved 
exactly Here we just quote the results for the trans- 



as already found in [1, Q. For the Duke exp., this be- 
havior is not dominant since it occurs only at very short 
times [smaller than 1 ms]. However, for times much 
longer than the trap period t ^ 1 there is a transi- 
tion to an energy dominated regime where (j3(x,p,t) ~ 
(j>{E, t) ~ il.oEt, where ftoE emerges as an effective 
transverse precession frequency, (ii) A second regime 
occurs for 1 <C t <C l/|f^o| where the phase difference 
(t){E,t) - (l)iE',t) ~ flotiE - i;') < 1 is smah - this is 
the regime of correlation between motional energy and 
transverse spin. Averaging over fast axial oscillations of 
small amplitude, we obtain the approximate result: 



m||(:r,i) « ^egn{x)'/'il - ^x')0=not' 



(10) 



This gives the dominant contribution to the spin segre- 
gation in the Duke experiment. Note that there is a sign 
inversion as compared to regime (i) and that the scaling 
is now t"^ instead of [lB| . The typical time for trigger- 
ing the spin segregation is ttrig ~ l/\/|.gr2o|- The dotted 
line of Fig. [2 shows a plot of (10) for t - l/|f^o|- (iii) 
The third regime is an ergodic regime: when t ^ 1/\^Iq\, 
the phase difference between high and low energy atoms 
becomes large t) - (t){E',t) ~ QotiE - E') > 1, 
therefore saturation occurs. However, the saturation is 
not captured by the linearized equations as it is actually 
related to norm conservation and to the non-linearity of 
the kinetic equations V^. By comparing the numerics 



to the analytical solution of the linearized equations, we 
find that they agree for the first ~ 100 ms. 

A simple model captures the essence of the phe- 
nomenon. Assume that we can divide the atoms into two 
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classes, the "hot atoms" with some motional energy (for 
instance SksT) and the "cold atoms" with some lower en- 
ergy (for instance kBT/2). All the atoms oscillate quickly 
in the trap so that the precession rate of their spins de- 
pends on the average of the longitudinal magnetic field 
along their trajectory [16]. Since the hot atoms go fur- 
ther away from the center of the trap, their transverse 
spin component rotates at a different rate than for the 
cold atoms. Now, when a hot atom crosses the cloud of 
cold atoms, which has a different spin direction, its spin 
experiences a rotation due to the ISRE in forward scat- 
tering. This gives to the hot atom a longitudinal spin 
component that depends on the sign of the ISRE, while 
the cold atoms acquire the opposite component. The ef- 
fect is cumulative: each time the hot atom crosses the 
cloud forwards and backwards, the rotations of its spin 
are additive. After some time, a large fraction of the hot 
atoms is transferred towards one spin state, a large frac- 
tion of the cold atoms to the other. A situation where the 
two classes of atoms have antiparallel longitudinal spins is 
stable as long as lateral collisions are ignored, since ISRE 
does not affect antiparallel spin directions, and since any 
phase space distribution that depends only on the energy 
is time invariant. Finally, because the atoms explore a 
region of space that depends on their energy, the spin 
state of the cold atoms is dominant at the center of the 
trap, the spin state of the hot atoms at the edges. 

To illustrate how this spatial separation arises from 



spin-energy correlations, we arbitrarily separate the 
atoms in two equal groups, with axial motional energy 
larger or smaller than ksT \n2. The difference in density 
profile between cold n< (x) and hot n> (x) atoms is then 
easily obtained in terms of the error function: 

n< - n> = n{x)[2Re{ETf[y/ln2 - x^/2]} - 1] (11) 

If the correlation c between longitudinal spin and mo- 
tional energy is not 100%, but say c = 60%, this result 
is reduced by the corresponding factor. The dashed line 
in fig. [2] is a plot of 0.6(n< — n>) according to (11); a 
comparison with the other curves shows that this model 
already gives a reasonable understanding of the profile. 

In conclusion, two major features of the new phe- 
nomenon are: a cumulative effect of the ISRE produc- 
ing a strong correlation between motional energies and 
longitudinal spin directions, and its stability over long 
times. The nature of the effect is more ballistic than 
hydrodynamic. The shape of the profile depends on the 
properties of the trap, not of the interactions: any phase 
space distribution of the internal states that depends only 
on the energy remains invariant under time evolution in 
the trap. ISRE therefore plays a crucial role in creating 
the profile, but not in maintaining it. 
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